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We discuss the spectrum of the QCD Dirac operator both at zero and at nonzero 
baryon chemical potential. We show that, in the ergodic domain of QCD, the Dirac 
spectrum can be obtained from the replica limit of a Toda lattice equation. At 
zero chemical potential this method explains the factorization of known results into 
compact and noncompact integrals, and at nonzero chemical potential it allows us 
to derive the previously unknown microscopic spectral density. 



1. Introduction 

Because of the spontaneous breaking of chiral symmetry and confinement, 
QCD at low energy is a theory of weakly interacting Goldstone bosons. In 
the spontaneously broken phase, the QCD partition function is a nonan- 
alytic function of the quark mass with a chiral condensate that is discon- 
tinuous as the quark mass crosses the eigenvalue axis of the QCD Dirac 
operator. The strength of this discontinuity is proportional to the eigen- 
value density, a relation known as the Banks-Casher formula ^. 

A theory with spontaneously broken chiral symmetry that is much sim- 
pler than QCD is chiral Random Matrix Theory ^. This is a theory with 
the global flavor symmetries of QCD in which the matrix elements of the 
Dirac operator are replaced by random numbers. Although, this theory 
is zero dimensional, chiral symmetry is broken spontaneously in the limit 
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of infinitely large matrices, and the mass of the non-Goldstone modes di- 
verges in the limit N ^ oo, where N is the matrix size. Therefore, in the 
thermodynamic limit, iV — > cxd, chiral Random Matrix theory reduces to a 
theory of Goldstone bosons for which, in the limit of small quark masses, 
the Lagrangian is just the mass term of the chiral Lagrangian. This is the 
main reason why Random Matrix Theories have been so successful in this 
context. 

One of the questions we have been asking is whether we can identify a 
parameter domain where QCD and chiral Random Matrix Theory reduce 
to the same theory of Goldstone bosons. The affirmative answer to this 
question is that this is the case if the Compton wavelength of the Goldstone 
bosons is much larger than the linear size L of the box. This requires that 
the quark masses nif <C F'^ /{TiL"^), which is an unphysical domain of QCD, 
so that the kinetic term of the chiral Lagrangian can be ignored (S is 
the chiral condesate and F is the pion decay constant.) However, even for 
realistic quark masses we can identify a parameter domain where QCD and 
chiral Random Matrix Theory behave the same, namely the domain where 
eigenvalues of the Dirac operator A ^ which is known as the 

ergodic domain ^. The reason is that the generating function of the Dirac 
spectrum is a QCD-likc partition function with additional ghost quarks 
with mass A. The condition for the validity of the Random Matrix Theory 
description of the Dirac spectrum is then that the Compton wavelength of 
Goldstone bosons composed out of ghost quarks is much larger than the size 
of the box. Indeed, such behavior has been observed in numerous lattice 
QCD simulations (as discussed in detail elsewhere ^'''). 

At nonzero chemical potential the eigenvalues of the Dirac operator are 
scattered in the complex plane. It has been shown that the Dirac spectrum 
remains in the ergodic domain if the inverse chemical potential is much 
larger than the size of the box ^. In this domain the Dirac spectrum can be 
described by a chiral Random Matrix Model that has been extended with a 
chemical potential However, until recently, this random matrix model 
had only been solved at the mean field level ^. The standard methods to 
derive nonperturbative results such as the supersymmetric method and 
complex orthogonal polynomial methods 12,13,14 j-^g^yg yjQ^ been successful 
in this case. So far the supersymmetric method failed because of technical 
problems in calculating the graded integrals, while the method of complex 
orthogonal polynomials failed because of the absence of an eigenvalue rep- 
resentation at non-zero chemical potential. The mean field analysis of the 
random matrix model was performed using the replica trick which was 
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widely believed to only work for the derivation of perturbative results 
However, it was shown recently that if a family of partition functions has 
certain integrability properties it is possible to obtain exact nonperturba- 
tive results by means of the replica trick I'^^i*. We will show that the ergodic 
limit of the phase quenched QCD partition function at nonzero chemical 
potential has the required integrability properties so that the hierarchy 
of phase quenched partition functions with a different number of flavors 
are related by a recursion relation which is known as the Toda lattice 
equation. The spectral density is then obtained from the replica limit of 
this Toda lattice equation. 

In the first part of this lecture we discuss the Dirac spectrum at zero 
chemical potential. We show that in the ergodic domain the spectral density 
can be obtained from the replica limit of the Toda lattice equation. This 
result explains a factorization property of the resolvent. In the second half 
of this lecture we study the quenched Dirac spectrum at nonzero chemical 
potential. Using the replica limit of the Toda lattice equation we derive the 
analytical result for the microscopic spectral density in the ergodic domain. 



2. The Dirac Spectrum in QCD 

The eigenvalues {Afc} of the anti-hcrmitian Dirac operator are determined 
by the eigenvalue equation 

D(j)k ^ iXkfpk- (1) 

Because of the axial symmetry the nonzero eigenvalues occur in pairs ±Afe. 
The number of zero eigenvalues is almost always equal to the topologi- 
cal charge of the gauge field configuration. The average spectral density, 
defined by 

p{\) = (Y.5{\-\k)y (2) 
can be obtained from the discontinuity of the average resolvent 

where, in both cases, the average is over gauge field configurations dis- 
tributed according to the QCD action. This can be seen by considering the 
rectangular contour in Fig. 1. Assuming that the average spectral density 
does not vary substantially along this contour, the average total number of 
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z = 




Figure 1. A typical Dirac spectrum (left) and the average resolvent G{z) in units of 
SV compared with lattice QCD data. 



eigenvalues inside this contour is p{X)l, where A is a point on the imaginary 
axis inside this contour. Therefore, if we integrate the resolvent along this 
contour, we obtain 



G{z) = il{G{iX + e) - G{i\ - e)) = 2Tiip{\)l. 



Using the symmetry of the spectrum we find ^ 



Re[G(a + e)] = T:p{X). 



(4) 



(5) 



In QCD, the chiral condensate defined by S = G{€)/V where V is the 
volume of space time, is nonzero because of the spontaneous breaking of 
chiral symmetry. The average level spacing near zero is therefore given by 



P(0) 



(6) 



which is also the scale of the smallest nonzero eigenvalue. Asymptotically, 
for large A, the Dirac eigenvalues approach the spectrum of a free Dirac 
operator with eigenvalue density given by ~ . 
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Because the eigenvalues near zero are spaced as IfEV it makes sense to 
introduce the microscopic spectral density ^ 

In the ergodic domain of QCD this is a universal function that can also be 
derived from chiral Random Matrix Theory. 



2.1. Ergodic Domain of QCD 

The low-energy limit of QCD for Nf fcrmionic flavors is described by a 
theory of weakly interacting Goldstonc bosons parametrized by the unitary 
matrix U{Nf). If their mass 

(with F the pion decay constant) the kinetic term of the chiral Lagrangian 
factorizes from the partition function and the mass dependence of the par- 
tition function in the sector of topological charge v is given by ^'^ 

Z^^(M)^ f det-^C/e^^^TriA/cz-i+Mtc/]^ 

Ju{Nf) 

where the quark mass matrix is defined by M = diag(mi, • ■ • ,mN^). This 
is the ergodic domain of QCD. 

As we will explain below partition functions with bosonic quarks are es- 
sential for obtaining the resolvent and average spectral density of the Dirac 
operator. For bosonic quarks, the Goldstone bosons cannot be parameter- 
ized by a unitary matrix. The reason is that symmetry transformations 
have to be consistent with the convergence of the bosonic integrals. Let us 
consider the case of one bosonic flavor. Then 



dot M .™t *^ ) = -7 / c?^'/>irf^02exp 
la' m J TT^ 



(10) 



so that the exponent is purely imaginary for m = and convergent for 
Re(m) > 0. 

The most general flavor transformation of the action in (10) is a Gl{2) 
transformation that can be parameterized as 

U = e"V with H^=H and VV^ = 1. (11) 



February 1, 2008 19:3 Proceedings Trim Size: 9in x 6in 



minn06 



6 



For U to be a symmetry transformation for m ~ we require that 

-'Gr'o)-(.ro)' 

so that H has to be a multiple of 173 and V has to be a multiple of the 
identity. The V part of U is not broken by the mass term and is thus a 
vector symmetry. Only the symmetry transformation exp(s(T3) is broken 
by the mass term so that the axial transformations can be parameterized 

by 

f e" \ 

U=\q^-s] with sG(-cx),oo). (13) 
For Nf bosonic flavors the axial transformations are parameterized by 

A) ^^^^ ^'=^- 

The Goldstone manifold is thus Gl{Nf)/U{Nf). In the domain were the 
kinetic term of the chiral Lagrangian can be ignored, the effective bosonic 
partition function (indicated by the subscript ~Nf) in the sector of topo- 
logical charge v is given by 21,22,23 

Z^_j,{mf)=f det''C/e^^^^[^^"'+*^'^l. (15) 

JGl{Nf}/U{Nf} 

The resolvent can be obtained from a supersymmetric generating func- 
tion that contains one additional fermionic ghost quark and one additional 
bosonic ghost quark, 

/ dct{D + z) 



G{z) = d. 



(16) 



\Aet[D + z'), 

Therefore the low energy limit of this generating function contains addi- 
tional ghost Goldstone bosons and fermions with mass given by 2zS]/F^. 

For z <^ i^^/EL^ the z-dependence of this generating function is given by 
21,22 

^Vi-i(^^)=/ def't/e^s^siYlA/c/^+Mtc/]^ (^7) 

JGl{Nf + l\l) 

with M ~ diag(mi,-- - ,m]\[j , z, z') and Gl{Nf + 1|1) are super-matrices 
with a unitary U{Nf + 1) upper left block, a GI{1)/U{1) lower right block 
and Grassmann valued matrix elements elsewhere. The number of QCD 
Dirac eigenvalues that is described by this partition function is of the order 

1^ = 
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This number increases linearly in Nc for Nc —>■ oo which was recently con- 
firmed by lattice simulations 

An alternative to the supersymmetric method is to use the replica trick 
to calculate the resolvent. It comes in two different versions: the fcrmionic 
replica trick defined by 

and the bosonic replica trick defined by 

G(z)= hm^logZ!!^ (z). (20) 

N f^O — ivy 

If we take the replica limit of the fermionic (19) or bosonic (20) partition 
functions directly, we will obtain a result that differs from the supersym- 
metric calculation For almost two decades there were no methods to do 
reliable nonperturbative calculations with the replica trick. In the next sec- 
tion we will show that these problems with the replica trick can be avoided 
if the take the replica limit of the Toda lattice equation. 



2.2. Toda Lattice Equation 

We now consider bosonic and fermionic partition functions with all masses 
equal to z which only depend on the combination, cf. (9) and (15), 

X = zJ:V. (21) 

The unitary integral in the fermionic partition function (9) can be evalu- 
ated by decomposing U = ydiag(e'^''')y^^ and choosing the {0k} and V as 
new integration variables. By expanding the Jacobian of this transforma- 
tion, given by OA;<i I 6xp(z0fc) — exp(i9i)\^ , the different terms factorize into 
products of modified Bessel functions which can be combined again into a 
single determinant. The final result is given by ^^'^^ 

Z-^^ix) = det [/,+,_,(a^)]fe,z=i,.. ,7V, ■ (22) 

By using recursion relations for the Bessel hmctions, this result can be 
rewritten as 

^N,i^) - ^N,il,-i) [(x9.)^-+'/.(x)],_,^„ ..._^^_^. (23) 

Next we use the Sylvester identity ^7,28 -^j^j^j-^ jg valid for determinant 
of an arbitrary matrix A. It is given by 

CijCpq — CiqCpj = det{A)Cij-pq, (24) 
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where the Cy are cofactors of the matrix A and the Cij;pq are double 
cofactors. By applying this identity to the determinant in (23) for i = j = 
Nf — 1 and p = q = Nf, we easily derive the Toda lattice equation 29,30,3i 

(xdx) \ogZj^^{x) = 2NfX — ' . (25) 

This equation has also been derived as a consistency condition for QCD 
partition functions It also occurs in other application such as for example 
in self-dual Chern-Simons theory 

In the case of bosonic quarks (15) the positive definite matrix is di- 
agonalized as t/ = V^diag(e*'" )F''^. Choosing the {sk} and V as new in- 
tegration variables, we obtain after expanding the Jacobian given by 
nfc<i(exp(sfc) - exp(s/))(exp(-Sfc) - exp(-s,)), 

Z'^j.^ix) = det [K,+u-i{x%^i^^ ...^^^ . (26) 

As in the fermionic case, this result can be written in the form of a r 
function 

Z-nM) = ^^vTT^det [{xdxf'''K^{x)]^^^^^ ,„^^^_^ , (27) 

where we have used that and {—lYK^, satisfy the same recursion re- 
lations. Therefore, the bosonic partition function also satisfies the Toda 
lattice equation (25). 

The two semi-infinite hierarchies are connected by 

Jim^^(xa,)2logZ]^^(x). (28) 

By extending the Toda lattice hierarchy to include an additional spectator 
boson, it can be shown that "^'^ 
1 

m - 

N, 



lim —ixdx) log Z"^ (x) ^ \im xdxixdx +ydy) log Z'( _^{x\y) 



= xd^ lim 2;9j.logZi _i(x|?/) 

y^x 

= xdxxG{x). (29) 

Taking the replica limit of the Toda lattice equation (25) we thus obtain 
the identity 

xdxxG{x) = 2x^Z'({x)Z'i^{x), (30) 

which explains this factorization property. In the same way we can show the 
factorization of the susceptibility into a bosonic and a fermionic partition 
function 
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Inserting the expressions for Z\ and wc find 

G{x) = +i^,_i(x)/„+i(.T)) + -. (31) 

X 

For V = ^ this resuh is shown as G(x)IYjV by the sohd curve in the right 
figure of Fig. 1. Agreement with lattice data is found in the ergodic 
domain of QCD. 



3. Dirac Spectrum at Nonzero Chemical Potential 



At nonzero baryon chemical potential the Dirac operator is modified ac- 
cording to 



(32) 



This Dirac operator does not have any hermiticity properties and its eigen- 
values are scattered in the complex plane 36,37,38,39,40 _ p^j. gjjja,ll n we 
expect that the width of the cloud of eigenvalues ^ /i^. The average 
spectral density is given by 



p(A) = (^52(A-A.) 



(33) 



p{z,z*) = lim — dzdz' log Z„(2;, z*), 

n— »0 TTO 



and the average resolvent is defined as usual by (3). They are related by 

5..G(2)|,=A -^P(A). (34) 
The quenched spectral density is therefore given by the replica limit ^^'^^^^^ 

(35) 

with generating function given by (note that n counts pairs of quarks) 

Z„(z, z*) = (det"(D + + z)det"(-D + m7o + z*)) . (36) 

The low-energy limit of this generating function is a chiral Lagrangian 
which is determined by its global symmetries and transformation properties. 
By writing the product of the two determinants as ^ 



det(£) M7o + z)det(-i:» -I- ^70 + z*) 



id + II Q z 
id-Id z* 
z id'' + 











(37) 
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we observe that the U{2) xU (2) flavor symmetry is broken by the chemical 
potential term and the mass term. Invariance is recovered by transforming 
the mass term as in the case of zero chemical potential and the chemical 
potential term by a local gauge transformation. In the domain of fi and z 
where we can neglect the kinetic terms, the partition function is given by 



Z„{Z,Z*)^ f dt/e-^Tr[C/,S][C/-\S] + I^TrAf(C/+C/-)^ (33) 

Ju{2n) 



S = ( ^" M , f^^J" ° ) . (39) 



where 

,^ - . 

-In J V Z*l 

3.1. Integration Formula 

We have proved the following integration formula 

d[/det'';7e5T^[^^(^+^"')l+Sp°j'Tr[(^^'^"'^)''l 

U{2n) 

= det[(x9.0^(ya,)'zr(:.,y)]„<,,<„_,, (40) 

where 

/■I 

Z'[{x,y)^ AdA/,(Ax)/,(-A2/)e22^p"''=°''(2pcos-iA)^ (4^) 
Jo 

and c„ is an n-dependent constant. For example, consider the case n = 1 
and all Op = 0. Then the integral is given by Zi{x,y) which is a known 
integral given by 

^r(^,y)la,=o = f M\I,{\x)I,{~\y) 
Jo 

^ ylv{,x)l^-i{-y) + xly-i{x)l„{-y) ^^2) 

J.2 _ y2 ' 

which is a known result for the QCD partition function with two different 

44 

masses . 

3.2. Toda Lattice Equation 

Using the integration formula of previous section we find that the low- 
energy limit of the phase quenched QCD partition function at nonzero 
chemical potential is given by 

Zniz.z*) = ^-^^-^det [(.5.)nz*a..)'^r(-,^*)]o<M<n-i ' (43) 
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where 





Using the Sylvester identity as at zero chemical potential we obtain the 
Toda lattice equation 

zd,z d,. logZ„(z,z ) = —{zz ) ^^^^^ ^,^j2 • (45) 

The spectral density (35) follows from the replica limit of this equation. 
Using Zq[z, z*) = 1 we find the simple expression 

1 zz* 
p{z,z*) = hm — a,9,.logZ,':(z,2*) = —ZUz,z*)Z1,{z,z*). (46) 
n-+o nn 2 

What remains to be done is to calculate the partition function with one 

bosonic and one conjugate bosonic quark which will be completed in the 

next subsection. 



3.3. The Bosonic Partition Function 

In this subsection wc evaluate the low-energy limit of the QCD partition 
function at nonzero chemical potential for one bosonic quark and one con- 
jugate bosonic quark. Because of convergence requirements it is more com- 
plicated to derive the chiral Lagrangian in this case. By a careful analysis 
we find 

Z^_,iz,z*)= f dC/det''C/e-^^[^^^l[c/-,B]+^TrC-(c/-/c/-/)^ 

(47) 

where B is the baryon number matrix defined in (39) and the mass matrix 
C and the anti-symmetric matrix / are defined as 

Although this integral is convergent for Ime > 0, it diverges logarithmically 
for e ^ 0. We have also derived the partition function (47) starting from 
a chiral Random Matrix Theory at nonzero chemical potential and using 
the Ingham-Siegel integral 

The integral (47) can be evaluated analytically by using an explicit 
parameterization of positive definite 2x2 matrices. We find 

Z1,{z + ly) = C_i log e e ^.^-^ ( 4^2^,2 ) ■ (49) 
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The final result for the quenched spectral density is given by 



19 



XdXe~^^^'^'^'\I,{X{x + ty)VJ:)\\ (50) 



The constant has been chosen such that the /i — + limit of p{x, y) for large 
y is given by Y.V/tt (see below). The solid curve in Fig. 2 shows a graph of 
this result for y = and ^i^F^V = 16 in terms of the ratio 

versus xYiV ai y = Q. The dotted curve shows the result which is obtained 
when the Bessel function K^, is replaced by its asymptotic expansion. This 
result that was obtained from a nonhermitian eigenvalue model that is 
not in the universality class of QCD. An important difference between the 
two results is that the spectral density (50) for y = is quadratic in x for 
a; — > 0, whereas the result given by the dotted curve is linear in x for a; — > 0. 

Taking the thermodynamic limit at fixed z and the Bessel functions 
can be approximated by their asymptotic limit. This results in 

^^^'^^"i^^I^ I'-'l'^^^' ^^^^ 

and p{x, y) = outside this strip in agreement with a mean field analysis 
For the integrated eigenvalue density we then find 

dxp{x, y) = (53) 

I" 

in agreement with the eigenvalue density at /i = 0. 



4. Conclusions 

In this work we have analyzed the ergodic domain of the QCD partition 
function where the pion Compton wavelength is much larger than the size 
of the box. In this domain the QCD partition function reduces to a theory 
of weakly interacting Goldstone bosons for which the kinetic term in the 
chiral Lagrangian can be ignored. Independent of the quark masses, the 
generating function for the Dirac spectrum is in this domain for sufficiently 
small eigenvalues. 

We have shown that fermionic partition functions, bosonic partition 
functions and the supersymmetric partition function are connected by a 
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Figure 2. The quenehed spectral density at nonzero chemical potential in the ergodic 
domain of QCD (full curve). Also shown is a result derived from an eigenvalue model 
(dotted curve). The left hand plot is a zoom in of the right hand one. 



Toda lattice equation. This recursion relation makes it possible to derive 
nonperturbative results using the replica trick. In particular, this reveals 
the factorization of the resolvent and the susceptibility into products of 
simple bosonic and fermionic partition functions. In this article, we have 
considered the resolvent at zero chemical potential and the spectral density 
at nonzero chemical potential. 

The resolvent for quenched QCD at nonzero chemical potential ap- 
proaches zero linearly as a function of Rc(z). In the unqucnchcd theory, 
where the argument of the resolvent is also the mass in the fermion de- 
terminant, we expect a discontinuity in the resolvent as Iie{z) crosses the 
imaginary axis. Can we understand the differences between these two the- 
ories in terms of the spectrum of the Dirac operator? As is suggested by 
earlier Random Matrix Theory simulations '^^ there are significant differ- 
ences between the two. For example, in the unquenched theory, because of 
the phase of the fermion determinant, there is no reason that the spectral 
density is positive definite or even real. The first analytical results for the 
unquenched spectral density where recently obtained by James Osborn 
for a nonhcrmitian Random Matrix Model that is in the universality class of 
QCD at nonzero chemical potential. In the thermodynamic limit his results 
show strong oscillations but the connection with broken chiral symmetry is 
still a mystery. We hope to address this issue in a future publication 
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The replica limit of the Toda lattice equation is a powerful method that 
is also applicable to other partition functions with an integrable structure. 
For example, we mention the Ginibrc ensemble parametric correlations 
and the two-point function of the Gaussian Unitary Ensemble. Our expe- 
rience tells us that all universal results that can be derived from a complex 
random matrix theory (/3 = 2) can also be obtained from the replica limit 
of a Toda lattice equation. 
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